Regularity, asymptotic behavior and partial 
uniqueness for Smoluchowski's coagulation 

equation 



Stephane Mischler * 
CEREMADE, Univ. Paris-Dauphine 
Place du Marechal de Lattre de Tassigny 
75775 Paris CEDEX 16, France 

Jose Alfredo Canizo t 
Departament de Matematiques 
Universitat Autonoma de Barcelona 
08193 Bellaterra (Barcelona), Spain 

May 2009 



Abstract 

We consider Smoluchowski's equation with a homogeneous kernel 
of the form a(x,y) = x a y@ + x@y a with —1 < a < (3 < 1 and A := 
a + (3 £ (—1,1). We first show that self-similar solutions of this 
equation are infinitely differentiable and prove sharp results on the 
behavior of self-similar profiles at y = in the case a < 0. We 
also give some partial uniqueness results for self-similar profiles: in 
the case a = we prove that two profiles with the same mass and 
moment of order A are necessarily equal, while in the case a < we 
prove that two profiles with the same moments of order a and j3, and 
which are asymptotic at y = 0, are equal. Our methods include a 
new representation of the coagulation operator, and estimates of its 
regularity using derivatives of fractional order. 
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1 Introduction 

1.1 Smoluchowski's equation 

Smoluchowski's coagulation equation is a well-known model for irre- 
versible aggregation processes involving a set of particles which can 
join to form groups of two or more of them, which we call clusters. 
Below we will briefly present the model, but we refer the reader to the 
reviews [H O [10J [11] for a more detailed background on the equation. 

We study this equation in a continuous setting, meaning that the 
size of a cluster may be any positive number y € (0,oo). We are 
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interested in the time evolution of the density of clusters of each pos- 
sible size, given by a function / = f(t,y) which depends on the time 
t and cluster size y. The mass of / at a given time t is given by its 
first moment J* °° y f(t,y) dy. The continuous Smoluchowski equation 
reads: 

d t f(t,y) = C(f(t),f(t))(y), (1) 

where 

C(f, f)(y) - = l f a(z, y - z)f(z)f(y - z) dz 

(2) 

a(z,y)f(z)f(y)dz, 

where the coagulation kernel a = a(x,y) is a given nonnegative sym- 
metric function defined on (0, +oo) x (0, +oo), which governs the time 
rate at which a cluster of size x and cluster of size y aggregate. We 
write C(/, /) to emphasize that C is a quadratic operator, and to agree 
with later discussions in which we consider its associated symmetric 
bilinear form C(f,g). 

In many physical models the coagulation kernel a is a homogeneous 
function [Ij in the sense that for some A > it holds that 

a(hx, hy) = h x a(x, y) for x,y,h > 0. 

For our results we will take that as a hypothesis, and in fact we will 
assume that a is of the following form: 

a(x, y) := x a y^ + x^y a for x, y > 0, (3a) 
-Ka</3<1, X:=a + /3 G (-1,1). (3b) 

In some places, we allow a to be a finite linear combination of such 
terms. 

When the kernel is homogeneous, one may look for self-similar 
solutions, also called scaling solutions: solutions which are given by 
a rescaling of some fixed function g for all times. In other words, a 
solution / = f(t,y) to Smoluchowski's equation is self-similar if there 
exists some nonnegative function g such that 

f{t,y)=q(t)g(p{t)y) forallt,y>0 

for some functions q(t),p(t) > 0. Such a function g is called a self- 
similar profile (or scaling profile). In this work we will always consider 
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self-similar solutions with finite mass. If a is homogeneous of degree 
— 1 < A < 1, it is known (e.g., [1]) that for / to be a self-similar 
solution it must happen that 

q(t) = {t + t)~& and p(t) = (t + t)"^x 

for some constant to > 0, and that g satisfies the self-similar profile 
equation: 

2g + yd y g + (l-X)C(g,g)=0. (4) 

The existence of self-similar profiles has been established in [6, 4) for a 
general class of coefficients which in particular cover the case of a(x, y) 
satisfying ([3]). It is expected (and has been proved in some particular 
cases) that in general conditions solutions exhibit a universal self- 
similar behavior in the long time, meaning that they are eventually 
close to a self similar solution of the equation. The conjecture that this 
behavior does in fact hold is called the dynamical scaling hypothesis, 
and a lot of effort has been done to rigorously prove its validity. In 
particular, it is expected that with "regular" initial conditions (such 
as compactly supported data) solutions become eventually close to a 
self-similar solution with finite mass. For the special cases a(x, y) = 1 
or a(x, y) = x + y rather complete results have been proved [9], [T2j [13] . 
but for a general coefficient almost no information is available. It 
seems likely that further study of the scaling profiles is needed before 
a more complete understanding of dynamical scaling can be obtained, 
and our results aim in this direction. 

1.2 Description of the main results 

We first give a proof of the infinite differentiability of self-similar pro- 
files. This fact was already known in the case a < [1], and here we 
generalize this result to include the a > case: 

Theorem 1.1. Assume that the coefficient a is of the form ([3]), or 
is a finite linear combination of terms of that form, all with the same 
homogeneity degree A. Then, any self-similar profile with finite mass 
for Smoluchowski's equation is infinitely differ entiable on (0, oo). 

We remark that in the paper [7j , self-similar profiles for a = were 
shown to be C , and in [H Th. 4.3] solutions for a > were proved to 
have a C ' 9 Holder regularity for any < 9 < 1 — A. 
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Regarding the asymptotic behavior of scaling profiles at and oo, 
some estimates have been rigorously proved [H E] ) and are consistent 
with the very complete previous formal calculations by van Dongen 
and Ernst [T7] . Exponential decay of solutions as y — > oo was proved 
in [H [7] , but the rate of decay is still far from the rate expected by 
formal calculations in [T7j. As for the behavior at y = 0, one has to 
separately treat the a = 0, a > and a < cases: 

1. For a = 0, a very precise result was given by Fournier and Lau- 
rengot in [7J, making rigorous the conclusions in [T7J: a self- 
similar solution is asymptotic to a constant times y~ T as y — > 0, 
for some 1 < r < min{3/2,l + A}. An intriguing property of 
this result is that in order to know the numerical value of r one 
needs to compute the solution itself, as r is given in terms of the 
A-moment of the solution. 

2. For a > 0, the formal arguments in [17] suggest that a scaling 
profile g should be asymptotic to a constant times when 
y — > 0. So far, it has been rigorously proved that all the moments 
of order a > A of g are finite, while all moments of order a < A 
are infinite. 

3. For a < 0, exponential decay of solutions at y = was proved in 
[1]. Our contribution here is a refinement of this behavior which 
coincides with the formal result in [T7j, and which is needed in 
the proof of uniqueness for a < 0. 

Our precise result for a < is the following: 

Theorem 1.2. Assume that the coagulation kernel is of the form Q 
with a < 0. Then, if g is a (nonzero) self-similar profile for Smolu- 
chowski's coagulation equation, it holds that 

g(y) ~ K e~ A ^ (5) 

for some strictly positive constant Kq, where 

A(y) :=21o gy -(l-A) (^V + , 

and M a , Mp are the moments of order a and (3 of g. In addition, the 
function y i— > g(y)e A ^ is decreasing on (0, oo). 

Remark 1.3. Notice that in the case (3 > the term in y 13 inside A 
does not play any role in the behavior of g as y — > 0, 



5 



We do not prove the asymptotic behavior of self-similar profiles 
at y = in the case a = because this was already proved in [7J. 
However, the same techniques employed in this paper, in particular 
the rewriting of the equation in section [51 may be used to give a proof 
of it which is somewhat different from the one in [7]. Lemma [53] below 
easily implies the following additional information on a profile g: 

Theorem 1.4. Assume that the coagulation kernel is of the form © 
with a = 0. If g is a self-similar profile for Smoluchow ski's coagula- 
tion equation, then the function y i— > y T ~ 1 g(z) dz is decreasing on 
(0, oo), where 

r := 2 - (1 - X)M x [g}. 

We also prove the following partial uniqueness results: 

Theorem 1.5. Consider a coagulation kernel of the form ([3]) with 
a = 0. Assume that g\ and g2 are two self-similar profiles of the same 
mass for Smoluchowski's equation, and also that M\[gi] = M\[g2\- 
Then, g 1 = g 2 . 

Here and in the rest of the paper we use the notation M^[g] to 
denote the //-moment of a function g, for any [i € R: 

MM := / y ft g(y)dy. (6) 
Jo 

For coefficients with a < our result is: 

Theorem 1.6. Consider a coagulation kernel of the form ([3]) with 
a < 0. Assume that g\ and gi are two self-similar profiles for Smolu- 
chowski's equation such that 

M a [gi] = MM =: M a , (7) 
Mp[ 9l ] = M p {g2] =■ Mp, (8) 
lim gi (y)e A ^ = lim o 2 (y)e A ^, (9) 

where 

A(y) :=21ogy-(l-A) (^V + . 
Then, g 1 = g 2 . 

Remark 1.7. The limit at y = which appears in this result is proved 
to exist in section [6j see theorem 11.21 below. Using it it is easy to see 
that condition ([9]) in the above theorem is equivalent to the require- 
ment that g\ ~ g 2 when y — ► 0, this is, lim 2/ ^o5 l i(y)/52(y) = 1- 
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To our knowledge, no uniqueness result at all was available for self- 
similar profiles of Smoluchowski's equation. The natural conjecture is 
that the full uniqueness result should hold, this is, that two scaling 
profiles with the same mass are necessarily equal. This does hold in 
the better understood cases with a(x, y) = 1 and a(x, y) = x + y, but 
the techniques we use here do not seem enough to show this for general 
A. A central difficulty is the fact that the moments of the equation 
(M\ when a = 0, or M Q ,Mg when a < 0) are a global property of 
the solution that cannot be computed a priori (unlike the a(x, y) = 1 
or a(x, y) = x + y cases), and which are even not easy to approximate 
numerically [2]. 

1.3 Tools and method of proof 

Among the tools used to prove the above results we highlight a new 
representation of the coagulation operator, which is given in section 
[3l and a rewriting of the self-similar profile equation, given in section 
[H Let us comment on them briefly. 

In order to rewrite the coagulation operator, define a distribution 
associated to any function / of finite mass. We use the Banach space 
L\ of real measurable functions on (0, oo) with finite first moment: 

L\ := L l {{®, oo); y dy). (10) 

In general, we use the notation 

L\:= L\{Q,™);y k dy) (11) 

for k G R. 

Definition 1.8 (Distribution associated to / € L\). Given a function 
/ € L\, we define the distribution {/} on M. as 

/•oo 

({/}, 0) := / f(z) (<f>(z) - 0(0)) dz for G Cg°(K). (12) 
J o 

Note that when / is not absolutely integrable at 0, this is just the 
classical definition of the finite part of the integral J* °° 4>(y)f{y) dy [IB] . 
Here, we keep the same expression even when / is integrable. Then, 
for a coagulation kernel of the form ([3]), the coagulation operator may 
be written as 

C(g,g) = {y a g}*{y 9 g}. 
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Here we are considering C(g,g) denned as a distribution; for a more 
precise statement, see section [3l 

The above expression has the advantage of being simple and lend- 
ing itself to convenient and perfectly rigorous manipulations. For 
example, it is easy to recover the known expression for the primi- 
tive of C(g,g), as the convolution above commutes with integration 
and derivation operators. Furthermore, it gives some insight into the 
way C(g,g) works: if we are working, say, in the case a = 0, where 
g{y) ~ K y~ T when y — > for some K > 0, one can see that C(g,g) 
is in many respects analogous to a fractional derivative of g of order 
r — 1 , as a fractional derivative of this order is precisely the convolu- 
tion with {y~ T }, times a constant. This is crucially used in the proof 
of the differentiability of profiles in section [H where we carry out a 
bootstrap argument which shows that if a profile g is k times differ- 
entiable, then the self-similar profile equation implies that it must be 
in fact k + 1 — A times differentiable. The use of fractional derivatives 
comes naturally in this context, and actually a different version of this 
bootstrap argument has already been used in [7] to show that profiles 
for a = are C . There, it was necessary to first show a certain Holder 
regularity of the solution, and then use this information to obtain C 1 
regularity. Formulating this in terms of gain of fractional derivatives 
makes it easy to iterate the argument to obtain C°° regularity and 
extend it to other kernels. 

For the study of the behavior of solutions near y = and the proof 
of uniqueness we rewrite the self-similar profile equation by solving the 
differential part, along with any other term which can be separated 
and solved; see section [5] for a statement of this. 

Finally, let us sketch the idea of the proof of our uniqueness result. 
A fundamental obstacle that makes equation difficult to study is 
the fact that it involves the nonlocal term C(g,g), which makes it 
very different from an ordinary differential equation. However, the 
coagulation operator for a constant coefficient has a gain part which 
only uses values of the function g less than y, and a loss part which 
is nonlocal only through the appearance of the integral of g. The 
idea is then to solve the latter part of the equation, which is simpler, 
assuming the value of the involved moments is given, and then look 
at the remaining part as an equation which is local near 0, to which 
the kind of arguments used in the theory of o.d.e.s can be adapted. 

This idea works well for coefficients with a < 0, as then the so- 
lution decays rapidly near y = 0, but it is not directly applicable for 
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coefficients with a = 0, as then solutions are known not to be inte- 
grable near and the gain and loss parts cannot be separated in the 
same way (one cannot separate the integral of g, as this term does not 
make sense). But, one can still find a way to separate the equation for 
the primitive of g in a similar way, and then carry out the argument on 
it. This line of reasoning is followed in section [71 where theorems 11.51 
and ll.6l are proved. This idea depends crucially on the fact that in the 
case a < the self-similar solution behaves better than near 
y = 0, so that the linear operator L g (f) := C(f,g) is regularizing near 
y = 0. Hence, this method does not give new information in the case 
a > 0, where solutions are expected to be asymptotic to y _1_A [TT| . 

1.4 Organization of the paper 

In the next section we give the basic definitions of self-similar pro- 
file, and precisely define the coagulation operator as a distribution on 
(0, oo). In section [3] we give the representation of this operator in the 
way mentioned above, which in particular extends it naturally to a 
distribution on (0, oo). This result is a central idea in the proof of in- 
finite differentiability of scaling profiles, to which section |4] is devoted. 
Section [5] proves that the self-similar profile equation can be rewritten, 
as we briefly explained before, in a way which is very useful to prove 
the asymptotic behavior of a < profiles at y = (given in section [6]). 
Finally, we prove our uniqueness theorems 11.51 and 11.61 in section [7J 

In an appendix (section [SJ we include a very brief introduction to 
fractional derivatives in order to clarify the notation and definitions 
we are using, as they are not completely standard in the literature. 
We also prove some simple but delicate results on fractional differen- 
tiation which are used in this paper, and for which we could not find 
a reference which gives the explicit statement. 

2 Preliminaries: self-similar profiles 

When one wants to define precisely the concept of self-similar profile, 
it is a well-known inconvenience that in order for C(g,g) to be well 
defined by expression ([2]), g must meet certain conditions; the ones 
which are usually imposed are finiteness conditions on certain mo- 
ments near and oo, which are not satisfied by the natural solutions 
of (jij) , known to have a nonintegrable singularity at y = [7J H3 HI E] . 
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Hence, the definition of C(g,g) is often changed for a less restrictive 
weak formulation by integrating against a suitably regular test func- 
tion (ft with compact support on (0, oo), which we will do next. Also, 
C(g,g) is quadratic in g, and we will need to consider its associated 
symmetric bilinear operator, so we actually give a weak definition for 
the latter: 

Definition 2.1. Let a be a symmetric nonnegative measurable func- 
tion defined on (0, oo) x (0, oo), and /, g be locally integrable functions 
defined on (0, oo) such that 

/ / a{y,z)y z\f(y)g(z)\ dydz < oo. (13) 
Jo Jo 

We define the coagulation operator C(f,g) associated to the coagula- 
tion coefficient a as the following distribution on (0, oo): 

(C(f,g),<j>) 

oo ;»oo 

/ a(y,z)f(y)g(z)((j)(y + z) - (ft{y) - <f)(z))dxdz 
Jo 

for (ft € C^°(0, oo). (14) 

Remark 2.2. For (ft € C^°(0,oo) there is always a constant M > 
(which depends on the compact support of (ft) such that 

\<P(y + z ) ~ <i>(y) - fit 2 ) \<Cyz, for y, z > 0, 

which can be readily seen by the Mean Value Theorem. Then, the 
integral (fT2|) is finite by (fT^ , 

When / is regular enough, (C(f, f),(ft) can be see to be equal 
to Jo°° < A(y) C '(/'/)(y) d y' where C(f,f)(y) is given by expression flS}. 
Condition (|13f) can be somewhat weakened in some cases by loosening 
the integrability condition on / at +oo, but the above one is simpler 
and will be enough for our purposes. 

Definition 2.3 (Self-similar profile). Assume that the coagulation 
coefficient a is homogeneous of degree A. A nonnegative locally inte- 
grable function g : (0, oo) — > [0, oo) for which 

roo 

I yg(y)dy<oo, 
Jo 

poo poo 

/ / a(y,z)y z\g(y)g(z)\ dydz < oo 
Jo Jo 
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is a self-similar profile for Smoluchowski's coagulation equation if 
equation holds in the sense of distributions on (0, oo); this is, if 



2g + yd y g + (l-X)C(g,g) = 0. (15) 

3 A representation of the coagulation 
operator 

In this section we want to give a representation of the coagulation op- 
erator when a has the form ([3]) which sheds some light on its structure, 
and in particular will be very helpful to prove our regularity results 
later. To begin with, let us give a natural extension of C(f,g) from 
definition 12, II to a distribution on R: 

Definition 3.1. Take a coagulation coefficient a and functions f,g 
in the conditions of definition 12.11 We define the coagulation opera- 
tor associated to the coagulation coefficient a, applied to /, g, as the 
distribution C(f,g) on R given by 

(C(f,g),<j>) 

OO POO 

/ a(x, z)f(x)g(z)(4>(x + z) — 4>(x) — 4>(z) + </>(0)) dx dz 
Jo 

for all (j) G Cg°(R). (16) 

It is easy to see that this is well defined as a distribution on R, 
even with the weak requirements on / and g, and that it extends 
definition ^. 11 Note the addition of the term ^(0), which does not make 
a difference when (j) has compact support on (0, oo); later we will see 
how this extension comes about naturally. To give our representation 
for C(f, g) we will use the notation from definition [L8j Let us initially 
treat the case of a constant coefficient a = 1, and then see how it 
extends to other coefficients: 

Proposition 3.2. Take f,g € L\, and let Cq be the coagulation op- 
erator with a constant coefficient a = 1 as given in definition \ 3.1l 
Then, 

Co(f,g) = 2^J" * as distributions on M. 

Here, the convolution {/} * {g} is understood as a convolution of 
distributions with compact support to the left [16] . This expression is 
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surprisingly simple, and is well-suited for the study of the coagulation 
operator when the functions /, g have a singularity at 0, which is the 
case with some self-similar profiles. Its proof just consists of writing 
out the definitions: 

Proof. For <f> G Cg°(M), 

<{/}* {9} A) = <{/}, {n{g}) *cf>) , (17) 

(where 1Z is the reflection operator, 7Z(j)(y) := 4>(—y), defined by du- 
ality on distributions) and 



((K{g}) * 4>){x) = ({g}, r- x <f>) = / g{z){<j>{x + z) - <j>{x)) dz, 

Jo 

so from (|l?p we have 

POO 

= / /(^)((("^M) * 4>){x) - {{Tl{g}) * 0)(O)) 


/(x)( f°° g{z){d>(x + z)-<f,(x))dz 



g{z){<l>{z)-<l>{SS))dzj dx 

f(x)g(z) (<f>(x + z) - 4>(x) - <p(z) + 0(0)) dxdz, 

which is our result, in view of expression (j!6|) . □ 

This directly gives a representation of the coagulation operator 
C(f,g) with a coagulation kernel a satisfying ([3]), as the following 
relation holds for any /, g satisfying (113j) : 



C(f,g) = C (y a f,y f3 g) + C (y^f,y a g), 

where C(f,g) is the operator associated to a, and the equality is an 
equality of distributions on R. 

Theorem 3.3. Assume that the coagulation operator a is of the form 
(0), and take /, g which satisfy ([15]) . Then, the coagulation operator 
associated to a (as given in definitional^ can be written as 

C(f,g) = \ ({y a f} * {/<?} + {//} * {y a g} 

where equality holds as distributions on R. 
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We emphasize that this operator is defined as a distribution on R 
and, as pointed out after definition 13. 1\ extends the usual operator C 
from definition 12.11 (which is a distribution on (0, oo)). 



4 Local regularity 

With the representation of the coagulation operator C(g,g) given in 
theorem 13.31 the study of its regularity can be viewed as the study 
of the regularity of convolutions of the above type, which is a more 
manageable problem. The difficulty here is that when a > a scaling 
solution g is not integrable near 0, but is expected to be very regular 
locally, and hence we need to study the convolution of functions which 
have singularities at 0. Precisely, the following general result on the 
local integrability of scaling profiles near y = is known (see [6] or [3] 
for a proof): 

Proposition 4.1. Assume that the coagulation coefficient a is of the 
form ([3]), or is a finite linear combination of terms of that form, all 
with the same homogeneity degree A. Then, all self-similar profiles 
g for Smoluchowski's equation (in the sense of definition \2.3\) satisfy 
that 



For convenience, we will measure the regularity of a function by 
looking at how many of its derivatives are locally integrable on (0, oo). 
To study the regularity of C we will need to use an interesting relation- 
ship between the kind of singularity of a function near and the local 
integrability of its fractional integrals, which we give in lemma S3] be- 
low. Let us start with the following elementary lemma that we state 
without proof: 

Lemma 4.2. For < k < 1, 




y k g(y) dy < oo for all R > and all k > A. 




for z > 0. 



(18) 



Lemma 4.3. If f G L\, with < k < 1 



then 



D~ k {f} 



2 

< 




Li(M) ~ T(A; + 1) 



L 1 (0,oo) 
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Remark 4.4. We recall that we denote L\ := L 1 ((0, oo); y k dy) as in 
eq. CLE]). 



Proof. For any <fi £ C^°(M) we will prove that 

2 



D-«{f},< 



< 



T(k + 1) 



L!(0,oo) 



(19) 



which is equivalent to our inequality. We have: 

l-OO 

D- k {f}, 0) = ({/},£>_*$ = J f(z) (D„ k 4>{z) - £)_ fc 0(O)) cfe, 



(20) 



and the part inside the parentheses is 



D- fc <^) - D_ fc <K0) 
1 



r(fc) 



l Z" 00 

x)(x - z) fc_1 (fx — - / 4>(x)x k ^ 1 dx 



1 



T{k) 



(x) ((x - z) 



k - x - x^ 1 } dx 



m j 



x)x k 1 dx. 



(21) 



We put back this these two terms in (|20p , and bound them separately. 
The first term is when k = 1, and for < k < 1 we have 



fc-l Jfc-1 



|/(*)| / |0(x)| (x-^-x 



dx dz 



< 



\m\ (x-z) 



i 



fc-1 _ x k-l 



dx dz 



k 



z k \f(z)\dz, (22) 



thanks to lemma 14.21 As for the second term in (12ip , putting it into 
(1201) we have 



\f(z)\ I \</>(x)\x k - 1 dxdz < 



f(z)\ I x k - x dxdz 



z k \f{z)\dz. (23) 



Then, eqs. (p^ - (j23|) prove that 



D~ k {f}, 



< 



kT{k) 
14 



z k \f(z)\ dz, 



which proves inequality (|19p . taking into account that kT(k) = T(k + 
1). □ 

In the light of the above lemma, our next result can be understood 
as saying: if two functions are locally regular but have a nonintegrable 
singularity at 0, their convolution is slightly less regular. How much 
less regular it is depends on the nature of the singularity. In particular, 
the local regularity of the convolution depends only on local properties 
of the initial functions, which is a general property of the convolution 
operation. In the next lemma, the reader can keep in mind that \jl will 
be negative when we use it, so a function / for which D^f is integrable 
may well be not integrable near 0, as lemma [4731 makes clear. 

Lemma 4.5. Let T , S be two distributions on R with support on 
[0, oo) (this is, T, S £ L)' L ), and assume that 

1. For some D V T and D V S are locally integrable on (0, oo). 

2. For some n < v, D^T, D^S are locally integrable on R. 
Then, the distribution D^ +U (T * S) is locally integrable on (0, oo). 

Proof. We break T and S into a part near 0, an intermediate part, 
and a part near oo. For this, choose < e < 1/4. We can find smooth 
nonnegative cutoff functions <I>0; &2 on (0, oo) such that 

$o = 1 on (0, e), $ = on (2e, oo) 

$i = 1 on (2e, -), $i = on (0,e) U (-,oo) 
e e 

2 1 
$2 = 1 on (-, oo), $ 2 = on (0, -) 

and such that 

$o + $i + <h = 1 on (0, +oo). 

In other words, ^2 form a partition of unity on (0, oo) subor- 

dinated to the open cover (0, 2e) U (e, 2/e) U (1/e, oo). Then, 

S = S<5> + + S$ 2 =: S + S x + S 2j 
T = T$ + r$i + T$ 2 =: To + T x + T 2 , 

where we have denoted Si := S&i, Ti := T$j for i = 0, 1, 2. We can 
break the convolution S * T by using this decomposition. Note that 
for i = 0, 1, 2, S" 2 *Tj is zero on (0, 1/e), and the same happens with S, 
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T interchanged; as we are only interested in studying the regularity of 
S * T on a bounded interval, we can disregard these terms and write 

S * T = S * T + S * Ti + Si * T + S t * T x on (0, 1/e). 

Similarly, the term So * To is zero on (4e, oo), so we have 

S*T = S *T 1 + S 1 *T + S 1 *Ti on(4e,l/e). (24) 

Then, we can write D^ +u of each of these terms by using theorem 18.81 



By the hypotheses of the lemma, we can see that all of the terms 
that take part in the convolutions on the right hand side are inte- 
grable functions, as the product by <£o or $1 does not change their 
local regularity properties (theorem I8.10p . Let us do the reasoning for 
D^S\: as D v S is integrable on (4e, 1/e) by hypothesis, we have that 
D U (S&\) = D V S\ is integrable on that interval by theorem 18. 10) then, 
thanks to lemma [8.111 D^S\ is also, as \i < v. The rest of the terms 
can be treated analogously, and are seen to be integrable without the 
help of lemma 18.111 

Then, all the terms on the right hand side of (125p -( j27l) are con- 
volutions of integrable functions, and hence are integrable, and (|24p 
then proves that D fJ,+u (S * T) is integrable on (4e, 1/e). □ 

Proposition 4.6. Take a coagulation coefficient a which is of the 
form ([3]) with a > 0, and C the coagulation operator associated to a 
(given by definition \3.1\) . Assume that g : (0, 00) — > R is such that 

1. y k g is locally integrable on [0, 00) for some 1 > k > /3, 

2. and D v g is locally integrable on (0, 00) for some v > 0. 
Then, D a ~ k+u C(g, g) is locally integrable on (0, 00). 

Remark 4.7. The hypothesis that a > is given for convenience, 
as we only need the result in that case; however, the lemma is true 
and proved in the same way also for negative a with the additional 
requirement that 1 + a > k > (3, so that {y a g} makes sense according 
to definition 11.81 



D^ +u (So * Ti) = {D^So) * (D v Tx) 
D^ L+u {S l * T ) = (p v S x ) * {D^To) 
D^ +u {Si * Ti) = {D^Si) * (D"Ti). 



(25) 
(26) 
(27) 
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Remark 4.8. The result works also when o is a linear combination of 
terms of the form ([3]) with the same A. In this case, if cti,/3j are the 
exponents in the i-th term in the linear combination, the proposition 
is true taking (5 := maxj{/3j} and a := mhijjaj}. 

Proof. It is enough to prove it for a coagulation coefficient of the 
form ([3]) with a > 0, as then we can apply the result to each term 
of the linear combination. For such a coefficient, the representation 
theorem 13.31 shows we can write C{g,g) as 

C(g,g) = {y a g}*{y f3 9}- 

As y i— > y k g(y) is locally integrable on [0, oo), we know that 

• y k ~ a y a g is locally integrable on [0, oo), and 

• y k ~^y^S is locally integrable on [0, oo), 
and hence from lemma 1*01 

• D a ~~ k {y a g} is locally integrable on R, and 

• D@~ k {y@ g} is locally integrable on R, so D a ~ k {y 13 g} is also (as 
(3 > a; see lemma IB. lip . 

In addition, both D u {y a g} and D u {y^g} are locally integrable on 
(0, oo), as {y a g}, {y^g} are equal to the functions y a g, y^g, respec- 
tively, on that set, and then theorem 18.101 applies there. Hence, we 
obtain our result as an application of lemma 14.51 with S := {y a g}, 
T := {yPg} and fj, := a — k □ 

Now we can finally prove theorem ll.lt 

Proof of theorem \l.l[ Take any A < k < 1, with A the homogeneity 
degree of a. Then, y k g is locally integrable on [0, oo) as recalled in 
proposition 14.11 (actually, we know it is integrable) . We will show the 
following: if, for some v > 0, D v g is locally integrable on (0, oo), then 
jju+i+a-kg j s a j go locally integrable on (0, oo). Asv + l + a — k>v, 
this implies that g is infinitely differentiable by a bootstrap argument 
starting with v = (g is locally integrable by definition). 

To show this, write the equation for a self-similar profile as 

2g{y) + yD 1 g + {l- \)C(g,g) = as distributions on (0, oo). (28) 

Assume that D v g is locally integrable on (0, oo). Then, by proposi- 
tion [4761 D u+a ~ k C(g, g) is locally integrable on (0, oo), and hence eq. 
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(f28|) shows that D v+a ~ k {yD l g) is locally integrable on (0,oo). By 
theorem 18,101 the same is true of D u+a ~ k D 1 g = D v+a ~ k+1 g, which 
proves our claim. Hence, g is infinitely differentiable. □ 



5 Rewriting the self- similar equation 

One of the techniques that we use in order to study the behavior of 
scaling profiles near y = and the uniqueness of self-similar solutions 
is a way of rewriting equation (|15|) in which we "solve" the differential 
part of the equation as far as possible. We introduce this method 
next: 

Lemma 5.1 (Solution of an o.d.e.). Let g be an absolutely continuous 
function, and let h, /i be locally integrable functions, all of them defined 
on (0, oo). // the following equation holds 

Kv)a{y) + yg'{y) = Kv) f or almost all y > o, (29) 

then g is given by 

g(y) = K(y)e- A ^ for all y > 0, 
where A and K are absolutely continuous functions which satisfy that 

A'(y) = for almost all y > 0, (30) 

y 

K'(y) = ^e A{y) h(y) for almost all y > 0. (31) 

Proof. We remark that one may find the expression of g by the method 
of variation of constants. To prove the result, a direct check shows 
that if we take functions A, K satisfying (|30p . (|3ip and define 

g(y) := K{y)e-^\ 

then equation (|29h holds with g instead of g. We may add a constant 
to K so that g(l) = g(l), for instance. Now, if we regard (|29p as an 
ordinary differential equation for g, then both g and g are solutions of 
it in the sense of Caratheodory, and then general uniqueness theorems 
(see, e.g., [8]) prove that g = g. □ 
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Lemma 5.2 (Primitive of C). Assume that the coagulation coefficient 
a is of the form ([3]) with a = 0, and take a function g € L\ n L\ (this 
is, with finite mass and finite moment of order X). Then the primitive 
ofC(g,g) can be written as 

D^C(g,g) = G*(y x g)-M x [g}G, (32) 

where G is the function given by 

/•oo 

G{y) := / g(z) dz for y > 0, G(y) := for y < 0. (33) 
Jy 

We remark that this equality is an equality of distributions on R. 

Remark 5.3. This is not a new result: one may arrive at the same 
expression by taking a characteristic function of the interval (y, oo) 
as the function eft in . The main point of the lemma is that it is 
rigorously proved as an equality of distributions. 

Proof. The representation of C(g,g) in theorem 13.31 gives 

C(g, g) = {<?} * b x g} = {g} * (y x g) - M x [ g ]{g}, (34) 

as g has finite A-moment. Hence, taking the primitive D_i (see eq. 
(fTT2j) ) we have 

D^C(g,g) = * (y x g) - M x [g]D^{g} 

= G*(y x g)-M x [g]G, 

where we have used the derivation rule for a convolution (cf. theo- 
rem [878]) and the fact that, as distributions on R, 

D-!{g} = G. 

□ 

Lemma 5.4. Assume that the coagulation coefficient a is of the form 
([3]) with a = 0. If g is a self-similar solution of Smoluchowski's equa- 
tion, it holds that 

G(y)=y 1 - T K(y), (35) 
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for some absolutely continuous function K, where G is given by (|33p . 
and where 

t := 2 - (1 - X)M x [g] (36) 

K'(y) = -y T - 2 h(y) for almost all y > 0, (37) 

h:=(l-X)(G*(y x gfj, (38) 

We note that for the convolution in A38\) it is assumed that g(y) is 
for y < 0. 

Remark 5.5. Note that the quantity r here is in agreement with that 
in P3E]. 

Proof. It is known that in the case a = the moment of order A of 
a solution g is finite [HE]. Hence, using ([32]) we can rewrite equation 
([11) by taking its primitive: 

G - yg + (1 - A) (g * (y A 5 ) - M a [ 5 ]g) = 0, (39) 

which holds for all y > 0. Equivalently, 

(1 - (1 - A)M A [ 5 ]) G - 2/5 + (1 - A) (G * (y x g)) = 0. (40) 

Rewrite this as 

(r-l)G-yg + h = 0, (41) 

where r and /i are given by (|36[) and (|38|) . Now, if we solve for G in 
eq. (|4ip by using lemma ISTTI with the independent term h, we obtain 
([35]) and ([37]) in our result. □ 

Lemma 5.6. Assume that the coagulation coefficient a is of the form 
([3]) with a < 0, and take a self-similar solution g of Smoluchowski's 
equation. Then it holds that 

g(y) = K(y)e- A ^ (42) 

for some absolutely continuous function K such that 

K\y) = -e K ^h{y) for almost all y > 0, (43) 
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with 

h:=-(l-\)(y a g)*(yf 3 g), (44) 
A(y) := 2 log y - (1 - A) + ) (45) 

when /? > 0, anci 

A(y) := (2 - (1 - A)M a ) log y - (1 - A)^% a (46) 

a 

w/ien /? = 0. 

Proof. Let g be a self-similar profile for Smoluchowski's coagulation 
equation with such a kernel a. It is known [1] that it is infinitely 
differentiable and has finite moments of all orders, and hence satisfies 
equation (JI|) in a strong way: 

2g + yd y g + (l-X)C(g,g) = 0, 

or, separating the gain and loss parts of C(g,g), 

2g + yd y g+(l-\)(y a g)*(yPg) 

- (1 - A) Mg y a g - (1 - A)M a /<? = 0, (47) 

where M a and Mp are the moments of order a and /?, respectively, of 
g. Now, if we apply lemma [57T1 to equation (|47j) with the independent 
term /i given by (|4"4"|) and fi := 2 — (1 — A)(Mg y a + M a yP), we obtain 
that 

g{y) = K{y)e-^\ (48) 
for some absolutely continuous functions K, A such that 

A'(y) = - (2 - (1 - A)(M^ y a + M a /)) a.e. y > 0, (49) 

K'(y) = -e A(y) h(y) a.e. y > 0. (50) 

y 

We may actually choose A as a particular primitive, as then the in- 
tegration constant for K can be adjusted so that (|48h is still true. 
Then, we can take A as in ()45|) (or (|46p when /3 = 0), and the result 
is proved. □ 



21 



6 Asymptotic behavior at y = for 
kernels with a < 

In this section we prove theorem 11.21 Assume that the coagulation 
coefficient a is of the form ([3]) with a < 0. Let g be a self-similar 
profile for Smoluchowski's coagulation equation with such a kernel a. 
Then from lemma loTUI we know it holds that 

g(y)=K(y)e-^\ (51) 

for some absolutely continuous functions K such that (|43p - (|44p hold, 
and with A given by (|45|) (or (|46p when /3 = 0). 

We will prove that X is bounded on (0, R) for any R > 0. Obvi- 
ously, X is bounded on any interval (5, R) with < 5 < R, as is clear 
from (|5ip . so the point is in proving that K is bounded on (0,5) for 
some 5 > 0. 

Then, take e, <5 > and define 

N e := sup g(y)$ e (y), 
ye(p,S) 

where 

fe A ^ if < y < e 

The dependence of N e on e is explicitly noted because we want to take 
the limit e — > 0; of course, iV e depends also on 5, but we do not write 
the dependence explicitly, as our intention is to fix S at some value, 
which has not been chosen yet. We can give a bound for g in terms of 
N e : 

g(y) < N^iy)- 1 for y 6(0,(5). (53) 



(54) 



On the other hand, from (I43p . 

"5 



*T(y) = K(5) - (1 - A) f -e^ z) h(z) dz for y > 0. 



Let us find a bound for /i using (|53p and our knowledge that g is a 
bounded function [1]: 

<?(y) < K for y > 0. (55) 

We have 

/i = -(l-A)(y a 5 )*(2A). ( 56 ) 
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To bound this for y E (0,5), take 5 small enough so that e A ^ is 
increasing on (0,5) (and so is <J> e (y) _1 ), and then 

(y a g)*(y p g) < N t K ^(y)- 1 ( y a )*(y p ) = N e K x d> e (y)-y +1 , (57) 

for y G (0,5). Then, from ([56]) . 

[%)| <iVe^ 2 1>e(y)-V +1 , (58) 

and continuing from (|54p . taking into account that e^ z '<& e {_z)~^ is 
decreasing on (0,5), 

r s 

\K(y)\ < K(5) + (1 - X)N e K 2 / z A e A(z) $ £ (z)- 1 dz 

Jy 

r s 

< K(5) + iV,^ e A{y) ^e^)' 1 / z x dz 

Jy 

< K(5) + N e K 4 e A(s/) $ e (2/) _1 5 A+1 for y G (0, 5). (59) 

Hence, multiplying by e~^ y ^ e (y) (which is always less than 1), taking 
into account that g<fr e = K e~ A< 3> £ , and taking the supremum over 
(0,3/), 

N e <K(5) + N t K 4 5 X+1 , (60) 

As A + 1 > 0, taking 5 small enough gives a bound for N t which is 
independent of e, and hence proves that K is bounded on (0,5). This 
in turn implies that it has a strictly positive limit at y = 0, as it is a 
nonincreasing function, which can be seen from eqs. (I43j) - (|44p . 



7 Partial uniqueness of scaling profiles 

Let us prove theorems 11.51 and 11.61 on the partial uniqueness of scaling 
profiles. We will prove the following proposition, which already con- 
tains theorem 11.61 and which will be seen to easily imply theorem 11.51 
(see section [73]) : 

Proposition 7.1. Assume that the coagulation coefficient a is of the 
form (J3|) with a < 0. Assume that g 1 and g 2 are two self-similar 
profiles for Smoluchowski 's equation, and also that 

1. in the case a = 0, M\[gi] = M\[g2\, and 

hm^- 1 G 1 (y) = hm^- 1 G 2 (y), 
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where r is given by (|36j) . and G , G 2 are the primitives of g 1 , 
g 2 based at +00, defined as in (|33p. 

2. in the case a < 0, M a [gi] = M a [g 2 ], Mp[g\] = Mp[g 2 ], and 

lim g 1 (y)e A ^ = lim g 2 {y)e k{ - y \ 
y^O y^o 

where A is given by (|45p . 
Then, g 1 = g 2 . 

We prove this proposition in two parts: in the next section we 
prove that the result holds near y = 0, and in section 17.21 we prove 
that the result is global. 

7.1 Local result 

We will start by proving proposition 17. II locally near y = 0; this is, we 
will show that under the same hypotheses there exists 5 > such that 
g 1 (y) =g 2 (y) for < y < 5. 

7.1.1 Proof for a = 

Step 1: Rewriting the equation for a difference of profiles. 

Assume that the coagulation coefficient a is given by ([3]) with a = 
(so j3 = A). If g is a self-similar profile, then from lemma 15741 

hold, so 

K(y) = K - f z T - 2 h(z) dz, (61) 
Jo 

with Kq given by 

K := lim K(y) = lim y^Giy), (62) 

y^O y-*0 

which is known to exist and be strictly positive [7]. Note that here we 
have integrated eq. ([37|) between and y, which can be done once we 
know K has a limit at y = 0. 
Gathering the above, we have 

y T - l G{y) = K - (1 - A) z^ 2 (g * (y A 5 )) (z) dz, (63) 

which is a remarkable equation in the sense that it is local near 0, and 
the term to the right is more regular near than that to the left, as 
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will be precised later. We also remark that in (|63p the parameters r 
and Kq depend on g. 

Now, let us obtain the corresponding equation for the difference 
of two self-similar profiles. Let g\, g 2 be two solutions of self-similar 
profiles, and assume that they satisfy the conditions in proposition l7.11 
Then, both g\ and g 2 satisfy (|63p with the same r and Kq, and we 
can take the difference to get 

y T ~ l AG(y) = £ z T ~ 2 ((AG) * (y x gi )) (z) dz 

z- 2 [G 2 *(y x Ag))(z)dz, (64) 

where 

Ag:=g x -g 2 , AG := G\ - G 2 . 
Now, to obtain local uniqueness near 0, take e > and write 

N = N(e):= sup y T ~ l \AG{y)\ . (65) 

Observe that this quantity is known to be bounded thanks to [7]. Let 
us prove from equation (|64p that, if we take e small enough, then iV 
must be 0. 



Step 2: Estimate for the first term. Constants independent 
of e will be denoted by K%, K 2 ... We will use the following bound, 
which holds for all solutions g of (HJ) [7] (and in particular for g\ and 
52): 

g{y) < Kiy~ T for y > 
for some constant K\ > 0, which implies that 



G(y) < K 2 y x - T for y > 0, 
for some other constant K 2 . Then, from (|66p . 



(AG)*(y x 9l )(z) 



< 



<KxN 



\AG(x)\ (x - z) x \gx(x - z)\ dx 
x^ T (x - z) x - T dx < K 3 Nz 2+x - 2t . 



(66) 
(67) 
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Hence, the first term on the right hand side of (|64p can be estimated 
by 



z- 2 ((AG) * (yV)) {z)dz 



<NK H I z x ~ T dz = N K A z 1+x - T . 

(69) 



Step 3: Estimate for the second term. For the second term 
in (|64f) we need to make AG appear instead of Ag. We use integration 
by parts to write: 



z T - 2 (G 2 * (y X Ag)) (z) dz 

' £ ( y G 2 *(y x Ag)) (z) dz 



= y 

+ (2-r) 



So ZT ~ 3 Jo * (yA A5) ) (X) dX dz - (70) 

Here, the boundary term at y = in the integration by parts vanishes, 
which is a consequence of the bound below in eq. (175D . which we will 
show next for the term G 2 * (y x Ag) appearing above. Write: 

£ (G 2 * (y x Ag)) (z) dz = D' 1 (g 2 * (y x Ag)) (z) dz 

= (G 2 *D-\y x Ag))(y). (71) 

Also, 

D-\y x Ag)(y) = —y x AG(y) + A F z A ^AG(z) dz, (72) 

J 



so 



(z)\ dz 



D-\y x Ag)(y) < y x AG(y) +A F z x ~ x |AG(; 

Jo 

< Ny^ T+x + NX ( z x - T dz < K 5 N y 1+x - T . (73) 
Jo 

Hence, from (TTTI) . 

£ (G 2 *(y A A 5 )) (z)dz 



{G 2 *D-\y x Ag)){y) 



<NK 2 K 5 F z l - T {y- z) 1+x - T dz<N K % y 
Jo 



3+A-2r 



• (74) 
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And finally, gathering fl70|) and ([73]). 

J\ T - 2 (G 2 *(y A A 5 )) (z)cfe 

< N KQy T - 2 y 3+x - 2r + N K 6 |2 - r| [ V z t - 3 z 3+x - 2t dz 

Jo 

<NK 7 y 1+x - T . (75) 

Step 4: Final estimate Now, from eq. taking the supre- 

mum on (0, e) and using ([69]) and ([75]) one has 

N(e) < N(e)K 8 e 1+x - T . (76) 

Note that 1 + A — r > [7J, which is crucial for this argument, and is 
a particular property of the coagulation kernel we are using. Hence, 
for e > small enough, we have that 

N = N(e) = 0, 

and hence that 

9i(y) = 92(y) for < y < e, 

which proves the equality in proposition 17.11 locally near y = in the 
case a = 0. 

7.1.2 Proof for a <0 

Assume again that a is of the form ([3]), now with a < 0. Take two 
self-similar solutions g\ , g2 in the conditions of proposition 17.11 Then 
lemma 15.61 applies to both g\ and Q2 with the same A, and following 
the same reasoning as in the a = case we have 

Ag(y) = AK(y)e~ A ^\ (77) 

with 

{AK)'(y) = -e A ^Ah(y) for almost all y > 0, (78) 

y 

Ah = -(l-A)(y a Ag)*(yPg 2 ) 
-(l-\)(y a g 2 )*(yPAg). 

Here we denote AA := A\ — A 2 for any function A, and K\, hi,K 2 , h 2 
are the functions K, h associated to g\ , g 2 as in lemma 15.61 
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Then, from the equality of the limits in ([9]) we see that AK(y) — > 
as y — » 0, so integrating 



AK(y) = [ V -e A{z) Ah(z)dz. (80) 
Jo ^ 

Let us find a bound for A/i using Ai'T. Take e > and call 

iV:= sup \AK(y)\. (81) 



Then, for the first term in ([79]), 

(y a Ag) * (/ ff2 )| < iV (y a e" A ^) * (/ 52 ) 

< NK e~ A(y) (/) * (y a ) < iVifi e- Afe) y A+1 , (82) 

where we have also used that g is bounded on (0, e); of course, we have 
much stronger information on its behavior at y = 0, which has been 
used through the fact that AK is bounded and has a limit at y = 0. 
For the second term in (1791) a similar calculation shows that 



(y a g 2 ) * (/A 5 ) | < NK 2 y x+1 e~ A ^. (83) 
Putting together ([82]) and ([55]) we obtain a bound for A/i: 

[A%)| < Ar# 3 y A+1 e- A(2/) , (84) 
and continuing from (|80|) . 

|Aif(?/)| < iViT 3 / z x dz = NK A y x+1 . (85) 
Jo 

Finally, taking the supremum on (0, e), 

N<NK 4 e x+1 , (86) 

and then (as A + 1 > 0) choosing e small enough proves that N = 0. 
Hence, 

9i(y) = 92(y) forye(0,e). (87) 
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7.2 Global result 



We now extend the local result in the previous section in order to 
finish the proof of proposition 17.11 After this it will also be easy to 
prove theorem 11.51 



End of proof of proposition 7.1, Take two self-similar profiles g\ and 
g 2 satisfying the hypotheses of the theorem. The results in previous 
sections show that g\ = g 2 on some interval (0, e), for some e > 0. 
Following a strategy usual in uniqueness theorems for ordinary differ- 
ential equations, we will show that whenever g\ = g 2 on an interval 
(0, yo), there is a 5 > such that g\ = g 2 on (0, yo + 5). Together with 
our local uniqueness result, a well-known argument then shows that 
gi = g 2 on (0, 00). 

So, assume that g\ = g 2 on (0, yo) f° r some yo > 0. Take 5 > (to 
be fixed later) and define 



N:= sup \Ag(y)\ 

= sup \Ag(y)\. (89) 

y&(0,y +S) 

Choosing 5 appropriately, we will prove that N = 0. Take the differ- 
ence of the self-similar equation for g\ and g 2 to get 

yAg = AG+(l-X)D_ 1 C(Ag,g 1 +g 2 ). (90) 
Hence, taking the supremum on (yo, yo + <^)j 

y N<\\AG\\ Lao(yojyo+s) 

+ (1 - A) \\D^C(Ag, gi +g 2 )\\ Lao{y0)yo+s) . (91) 

To bound the first term we have, for y € (yo, Ho + 5), 

ry 



|AG(y)| 



8g(z) dz 



.'/ii 



<N(y-y )<N(5-y ). (92) 



The bound for the second term in (|9ip depends on the type of kernel 
we are considering: 
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Bound for an a < kernel 

2C(Ag, 9l ) = {y a Ag} * + {y? Ag} * {y a 9l } 



(y a Ag) * {yP 9l } + {y? Ag) * {y a 9l }, (93) 



0, 



so, for 5 < j/ < 25, 
\C(Ag, 9l )\ 



< \\y a Ag\\ Lac ^ y) Mp + 



y"Ag 



L°°(0,y) 



M n < KqN 



where the first inequality has been obtained by applying Young's 
inequality for measures (or, alternatively, by writing out explicitly 
C(Ag,gi) with its classical formula and bounding each term). Hence, 



£\C(Ag, 9l )\ dz<K N(y-5). 



(94) 



Bound for an a = kernel 



2C(Ag, 9l ) = {y x Ag} * { 5l } + {Ag} * {y x 9l } 

= (y x Ag) * { 9l } + {Ag} * (y x 9l ) - M\{Ag}, (95) 

and then 

2\D^C(Ag, gi )\ 
< Dx^Ag) * D_ X { 91 } + £>_! {Ag} * (y x 9l ) + M A D_i{A0} 



D^^Ag) 



L°°(-oo,y) 

We need to bound ^Dx-ii^^g)^^^. 
1 /■» 



+ 2M A (y-<5)iV. (96) 



£> A _i(y A A 5 ) 



z A A 5 (z)(y-z)- A dz 



r(i - A) 



< K X N J\y - z)- x dz 



= K\N 



{y-5) l - x . (97) 



1-A 
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Finally, with (|94p and (|97p we can continue from (|9ip to get 

yoiV^KsiVd/-^) 1 ^, (98) 

so choosing 5 small enough we can deduce that N = 0. This finishes 
the proof. □ 



Let us finally prove theorem 11.51 

Proof of theorem \l.b\ Assume that a is of the form ([3]) with a = 0, 
and take two solutions g\, g 2 of equation m the conditions of 
theorem ll.5l this is, 



ygi(y)dy = / yg 2 (y)dy, 
o J o 

OO /"OO 

y X 9i(y)dy = / y x g 2 (y)dy. 
'o Jo 

Then, for any /i > 0, the function 

9i{y) ■= fJ- 1+X 9i(^y) for y > 

is another solution of equation @, which is a simple consequence of 
the homogeneity of the coagulation coefficient a (see, for example, \7\). 
We also check easily that 

00 /'OO 

y x g 1 (y)dy = / y x g x {y) dy (99) 

JO 

1 /*oo 

yg 1 {y)dy = ^ x - 1 ygi(y)dy (100) 
Jo 

/"OO /'CO 

lirn^- 1 / g 1 (z)dz= f , 1+x - T limy T ~ 1 gi(z)dz, (101) 

Jy 2^° J,, 

where r := 2 — (1 — \)M\[gi]. Hence, the moment of order A of g\ is 
the same no matter which \x we take, so we can choose \x > in such 
a way that 

/*oo /'OO 

/ y X 9i(y)dy = / y x g 2 (y)dy 
Jo Jo 

poo roo 

limy T_1 / g\{z)dz = limy T_1 / g 2 (z)dz. 
y^o J y y^o J y 

Hence, with this value of fi, by proposition 17.11 we have that g\ = g 2 , 
and in particular their masses are equal; as the masses of g\ and g 2 
are also equal, from (jlOOp we deduce that in fact [i must be equal to 
1, so g\ = gi = g 2 , which shows the result. □ 
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8 Appendix: fractional derivatives 



The extension of the concept of integration and differentiation to in- 
clude derivatives and integrals of noninteger order is a well established 
theory [16j 113 E! • Here we give a brief but self-contained introduc- 
tion to it and state without proof the main standard results. Some 
particular properties needed in the rest of this paper and which are 
not commonly encountered will be given with complete proofs below. 

For our purposes, the simplest and most general definition of frac- 
tional derivatives is given in the context of distributions, and can be 
found in the book by Schwartz [161 VI. 5]. The reader can check that 
our definitions are the same as those given there, even if the presen- 
tation is somewhat different. Other expositions are found in [151 114j . 
and we refer to those sources for the proof of the main results on 
fractional differentiation given below. 

In the following we will use the space C|j? consisting of all infinitely 
differentiable functions / : R — > R which have support bounded below; 
this is, those / which have support contained in [a, +oo) for some 
a € R. The notation is intended to suggest that the important part 
of a function / is to the right, if one represents the real line as usualQ 
Analogously, we define C£° as the set of all infinitely differentiable 
functions / : R — > R whose support is bounded above. 

We will first define fractional derivatives for smooth functions, 
and then extend the concept to distributions with a common dual- 
ity method. 

8.1 Fractional derivatives of smooth functions 

Definition 8.1 (Left fractional derivatives). For / S and real 
k > 0, we define the left integral or order k of f as 



where F is the Gamma function. For k = we just write D k f = f. 
For real k > we write k as k = n — s, with n > an integer and 

1 We prefer to write C£? instead of T> + , used in Schwartz's book, as in the literature 
related to Smoluchowski's equation the subscript '+' is frequently used to denote that 
functions in the involved space are nonnegative, which could cause confusion here, where 
we want to stress a property of their support. 




for y E R 



(102) 
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< s < 1, and define the left derivative of order k of f to be 

D k f:=—(D^f). (103) 

The above is a usual definition of fractional integrals and deriva- 
tives |15| . sometimes called the Riemann-Lebesgue definition. Names 
given to the above also differ slightly from place to place: we may refer 
to D k f as the left derivative of order k of f, or just the k-th derivative 
of f, for any real k (even for k < 0), thus emphasizing that all D k are 
part of a family of operators with common properties. 

There is a completely analogous concept of right derivative where 
integrals are taken from +oo: 

Definition 8.2 (Right fractional derivatives). For / G C£° and real 
k > 0, we define the right integral or order k of f as 

1 f'°° 

D - kf{y):= v(k) I f ( z K z - y } k ~ 1(iz for y GM > ( 104 ) 

where T is the Gamma function. For k = we just write D^f = f . 
For real k > we write k as k = n — s, with n > an integer and 
< s < 1 , and define the left derivative of order k of f to be 

d n 

Dkf ■= (-1)*^(£>-./). (105) 
Some easy consequences are the following: 

1. For k G R and / G Cg 3 , L> fc / is again on C|?, and the analogous 
result holds for and 

2. For integer k > 0, D fc is just the usual fc-th derivative of /, 
while D~ k is the fc-fold iteration of the primitive based at — oo; 
in particular, D~ l f is the only primitive of / which is at — oo. 
Analogously, is the usual k-th derivative of /, times (— l) k 
(see next for a natural reason for this definition), while is 
the fc-fold iteration of the primitive based at +oo. 

3. Right derivatives are the concept symmetric to that of left deriva- 
tives under the reflection of R: if we define the reflection of a 
function / : R — > R to be the function IZf : R — > R given by 

:= f(-y), then 

D k (f):=lZ{D k {lZf)) for fcGR,/GC£°. (106) 
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Alternatively, one can take this as a definition of right derivatives 
from the perhaps more natural concept of left derivatives. Note 
that the alternating sign in equation (|105p is unavoidable if we 
want to conserve this symmetry property. 

Remarkably, the following composition result holds: 
Theorem 8.3. For real k,j, 

D>(D k f) = D> +k f for any f € Cf , (107) 
Dj(D k f) = D j+k f foranyfEC?. (108) 

A proof follows from elementary analysis arguments. For this result 
to hold it is essential that our definitions 18.11 and 18.21 above have picked 
specific primitives (those which are at -co or +oo, respectively) out 
of all the possible primitives of a function /. Said in another way, 
the spaces C^?, C£° in which we are working only contain one of 
all the possible primitives of a given function /, and thus the above 
composition rule can hold. 

Also, it is easy to see that D k is the dual of D k in the following 
sense: 

Lemma 8.4. For f E and g G C£° it holds that 

/+oo r+oo 
D k f(y)g(y)dy= / f(y)D k g{y)dy. (109) 
-oo J — oo 

This suggests the definition for distributions given below. 



8.2 Fractional derivatives of distributions 

Consider the set V' L of distributions on R which have compact support 
bounded below. One can show that T>' L is the dual of C|° when the 
latter is equipped with a natural topology (which extends that of Cq°) 
[161 VI. 5]. So, T>' L should be thought of as (C^ )', which is useful for 
remembering that distributions in this space have support contained 
in (a, oo) for some a € R. For these distributions, one can define (T, ip) 
for any t/j € as 

(T,if>) := (T,$) (110) 

for any ip G Cq°(R) which coincides with tp on the support of T. Of 
course, this definition does not depend on the particular extension 
chosen. We define T>' L analogously, and also the pairing (T, cj>) for any 
T € T>' L , (j) E Cr. 
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Definition 8.5 (Fractional derivatives of distributions). Take k G R. 
For a distribution T G T>' L we define the distribution D k T as 

{D k T, ip) := (T,D k rP) for ^ G Cf (R). (Ill) 

Analogously, for a distribution T G we define the distribution D k T 

as 

(D k T,^) := (T, D k p) for G (R). (112) 

Here, D k <p and -DfcV' are the right and left fractional derivatives, re- 
spectively, defined in section 18.11 Note that the duality products here 
are well defined as indicated in (IllOj) . Also, this agrees with defini- 
tions [8J] and [8T2] when T is a function in C R (or C£°), as can be seen 
from (fTn9l) . 

Then, the D k are linear operators for which the composition rule 
(fTUTjl still holds: for any j, k G R, 

D j (D k T) = D j+k T foranyTeP^ (113) 
Dj (DkT) = D ]+k T for any T G T>' L . (114) 

The convolution of two distributions in D' R (or two distributions 
in T>' L ) is well defined [16, VI. 5], as we show below: 

Definition 8.6 (Convolution of a distribution and a smooth function). 
Given T G T>' R and <p G C|° we define T * <p as the distribution in Z?^ 
given by 

(T * cj), ip) := (T, (TZ(p) * ip) forVGCg 3 . (115) 

Note that (K(p) * ip £ C R . For T G £^ and the convolution 

T * ip is defined analogously. 

The convolution with a function in C R (or C£°) is regularizing, as it 
happens in the more familiar case of convolution with a C°° function 
with bounded support: if T G T>' L and ip G C^P, then one can prove 
that T * ip is equal to a function in C R , given by 

T*iP(y) = (T,T y iP) foryGR, (116) 

where (T y ip)(x) := ip(x — y) is the translation of ip by y. Similarly the 
convolution T * <p for T G T>' R and <p G C|° is a function in C|° . 
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Definition 8.7 (Convolution of two distributions). Given T,S G T>' R , 
we define T * S as the distribution in T>' R given by 

(T * S, tp) := (T, (KS) * %l>) for r/> G C£. (117) 

(As remarked before (|116p . ilZS) * V G The convolution T * S 

for T, S" G T>' L is defined analogously. 

The following well-known result on the derivation of a convolution 
holds in complete generality with these definitions: 

Theorem 8.8. For any T, S G 2?^ and any 

D k (T*S) = (D k T)*S. 

Although fractional derivation and integration operators of non- 
integer order are not local, their regularity properties nevertheless 
depend only on local properties of the function they act on, as can be 
easily proved by observing that D k f can be written as a convolution 
of / with a distribution which is C°° away from \16\ VI. 5, p. 174]. 
A manifestation of this is the following result: 

Theorem 8.9. Let T be a distribution on M with compact support to 
the left (this is, T G V'J, and such that T is C°°{U), for U C R a 
given open set. Then, D k T G C°°(U) for all k G R. 

It is well known that multiplying a distribution by a C°° function 
preserves its local regularity, and this is still true when one measures 
this regularity in terms of the integrability of a given fractional deriva- 
tive. As this result is not easily found in the literature and its proof 
is not obvious, we give it in section [8731 below: 

Theorem 8.10. Let T be a distribution on R with compact support 
to the left (this is, T G L)' L ), and assume that D^T is locally integrable 
for some /i > 0. Then, for any smooth function <i> on R, D^ l (^T) is 
locally integrable. 

Of course, the analogous result holds for a distribution with com- 
pact support to the right and right derivatives D^. 

Let us also prove a result which is used in this paper, which says 
that the property that D k T is locally integrable is stronger the higher 
k is: 
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Lemma 8.11. Let T E T>' L be a distribution on R with compact sup- 
port to the left, and assume that D k T is locally integrable on R for 
some k € R. Then, £) k ~ m T is locally integrable for all real m > 0. 

Proof. It is enough to prove it for k = 0, as then the general result 
is obtained by applying this particular case to the distribution D k T, 
taking into account the composition law (|113|) . 

Then, to prove it for k = 0, take T a locally integrable function 
on R with support contained on (i?, oo). Fix a compact interval [a, b] 
with b > R. For a test function <j) £ Cq°(R) with compact support 
contained on (a, b) and any m > we have 

(D~ m T, <(>)-- 

and hence 

\(D- m T^)\ < \\4 

= H\ 



which proves that D 

8.3 Regularity of a product by a smooth func- 
tion 

In this section we give the proof of theorem 18.101 on the regularity of 
fractional order of a product by a C°° function. We recall its statement: 

Theorem 8.12. Let T E T>' L be a distribution on R with compact 
support to the left, and assume that D^T is locally integrable for some 
\i > 0. Then, for any smooth function $ on R ; D^{^T) is locally 
integrable. 

The proof is broken into several lemmas. The following one some- 
times serves as a weaker substitute for the rule of differentiation of a 
product: 



1 



Tim) 



1 

oo 



T(y) f ' <t>{z){z-y) m - l dzd yi 



R 



h T{y)\ \\z-y) m - x dzdy 



r(m) J R 

I fb rb-y 



T(m) 



T(y)\ z m ~'dzdy 
Jo 



1 T is locally integrable on (a, b). □ 
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Lemma 8.13. Take < k < 1. If (f>,ip £ C°°(R) and have compact 
support to the right, then the following equality holds: 

D fc (0D_ fc ^)(y) = <KvMy) 

sm(nk) [°°^_}_ f X < p'^ x _ z f^ z _ y yk dzdx _ 



vr J v x-y 



Proof. We have D k (cf)D- k ij)) = DiD k -\{^D- k tj)). Let us calculate 
this: 

D k - X (4>D- k il>){y) 
1 



z) / tp(x)(x-z) k 1 {z-y) k dxdz 



T{k)T{\ - k) 

poo px 

= K X I ij){x) / <f>(z)(x- z) k - 1 {z-y)- k dzdx, (118) 



v 



where we have set K\ := l/(r(A;)r(l — k)) for short. By a change of 
variables u = {z — y)/(x — y), the inner integral can be written as 

(z)(x-z) k - 1 {z-y)- k dz= [ (j)(u{x-y) + y)(l-u) k - 1 u- k du, 
y Jo 

so we can calculate the derivative in y of (1118j) and obtain 



D k {<j>D_ k iP){y) = -Kx— j tf>(x) j ^{z){x-zf~ L {z-y)- k dzdx 



v J y 



d roc rl 

-Ki— / ip(x) / 4>(u(x - y) +y)(l - u) k ~ 1 u~ k dudx 



dy Jy 



o 



r(fc)r(i - k) 

poo p 1 

-Ki j ip(x) / 4>{u{x - y) + y)(l - u) k u~ k dudx 

Jy Jo 

^(y)<t>(y) - K x / il>(x) / (j)'{z){x-z) k {z-y)~ k dzdx, 



y 



x-y Jo 



where B(k, 1 — k) is the Beta function for the parameters (k, 1 — k). We 
have used the well-known relationship between the Beta and Gamma 
functions, and have undone our previous change of variables. This is 
the expression in the lemma; note that K\ is the constant that appears 
there. □ 



38 



Lemma 8.14. Take k > and $ G C£°(R). T/ien, /or a// ^> G Cg°(R) 
with compact support contained on a fixed interval (—oo,b) it holds 
that 

\D^D k ^){y)\ <(fe-y)||V|looll*lloo for ally <b. (119) 
Proof. For y > b, D_ k (&D k _iif))(y) is zero. For y < b we have 

|£>_ fc (&D fc _i^)(l/)| 

- r(fc)r(i - k) f 1*^1 ^ - f) fc_1 /" l#»0 1 (*-*)"* 



i 



< ; - ; 11-011 ||$|l / / (z-yY~Hx- z)~ K dzdx 

< r^)r(i"- fc fc) {b ~ y) M ~ m °° = {h ~ y) m °° m °» ' 

taking into account the relationship between the Gamma and Beta 
functions. □ 

Lemma 8.15. Take $ G C£°(R) and < k < 1. Fix a < b G R. 

Then, for all ip G Cq°(R) compact support contained on a fixed 

interval (—oo,b) it holds that 

\\D k ^D^)\\ LOO{ajb) KKIWW^, 

where K > is a constant that only depends on k, $ and the interval 
(o,6). 

Proof. For = the statement is trivial. For < k < 1, we use 
lemma [8. 131 to write, for y < b, 

D k {<S>D^){y) = 

sm ( 7rA; ) /" ^( X )_L_ [ X &( z )(x - z) k (z-y)~ k dzdx, 



where we have also taken into account that ip(x) is for x >b. Here, 
the first term has the straightforward bound ||$Vlloo — Halloo Halloo- 
As for the second one, 

1 f x 

ip(x) / &{z){x - z) k (z - y)~ k dzdx 



x-y 



IOO M MOO , X _y 



b 1 



< ||rU||<[.'L / — / (x-z) k {z-y)- k dzdx 

= K{b-y)U\\\W 



oo || Moo 
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where the constant K is 5(1 + k, 1 — k). Using this for a < y < b 
proves the statement. □ 

Lemma 8.16. Take $ £ C£°(R) and k < 1. Fix a < b e R. Then, 
for all ip G Cq°(]R) u>i£/i compact support contained on a fixed interval 
(—oo,b) it holds that 

\\D k ^D^)\\ Loo{ajb) <Jf|C, 

where K > is a constant that only depends on k, $ anrf i/ie interval 
(a, 6). 

Proof. Lemma [8.151 proves this when < A; < 1, and we can prove the 
general case inductively: if the result is valid for a given fc < 1, then 

D fc _i($Z>i_ fc ^) = I>fc-i($Di(Z>_ fc ^)) 

= ZVi^i^-fcVO - D k _ x {iD^){D^)) 

= D k ($D_ k i>) - D fc _i((Z)i*)(Z)_ fc ^)). 

For the first term we can use our induction hypothesis, and the second 
one can be bounded thanks to lemma (I119p . This shows the lemma. 

□ 

Proof of theorem \8.10[ Take R G K so that T has support contained 
in (R,oo). Fix a compact interval (a, 6) with b > R. For a function 
■0 € C^°(M) with compact support on (a, 6) we have 

(D"(*T),^) = {&T,Drf>) = (T,*Z> M V> 

= (D-^D^T, SD^) = (D»T, D_ M *D M ^> . 

Taking into account that D^T has support contained in (i?, oo) and 
that D-^QD^ijj is smooth and has support contained in (—00,6), we 
have 



K£>"(*T),V>| < ll^|| L1(Ri6) P-m^mVIIxo-^) 

<K||OT|| il(ii)6) 

thanks to lemma 18. 161 □ 
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